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REPRESENTATION OF TURBULENT SHEAR STRESS BY A PRODUCT 
OF MEAN VELOCITY DIFFERENCES 
by Willis H. Braun 
Lewis Research Center 

SUMMARY 

A simple argument leads to the proposal of a quadratic form in the mean velocity for 
the turbulent shear stress. The proposed form is expressed as the product of two ve- 
locity differences whose roots are the maximum velocity in the flow and a "cutoff” ve- 
locity below which the turbulent shear stress vanishes. Application to pipe and channel 
flows yields the centerline velocity as a function of pressure gradient, as well as the 
velocity profile. The flat-plate, boundary-layer problem is solved by a system of inte- 
gral equations to obtain the friction coefficient, the displacement thickness, and the 
momentum-loss thickness. Comparisons are made with experiment. 


INTRODUCTION 

Theories of turbulence are generally classified according to the number of partial 
differential equations that must be solved in conjunction with the Reynolds equation for 
the time-mean flow (Reynolds, ref. 1; Launder and Spalding, ref. 2). If one additional 
equation (e. g. , an equation for the turbulent kinetic energy) must be solved, the theory 
is called a one-equation model; if, beyond that, a differential equation for the charac- 
teristic length of the turbulent motions must be solved, the theory is called a two- 
equation model; and so forth. 

The move toward more elaborate treatments of turbulent flows arises from the 
failure of the zero-equation models, that is, those which model turbulent stresses with 
an algebraically specified relation to the mean-flow quantities, to prove both general 
and accurate. The zero-equation theories have been constructed almost exclusively 
upon the hypothesis of Boussinesq (ref. 2) that the stress should be modeled as a product 
of an eddy viscosity and the mean-velocity gradient in the manner of Newton’s law for 
viscous fluids. Various forms of this hypothesis have been successful in describing a 
wide variety of flows of practical importance. Nevertheless, it is now generally agreed 



(Townsend, ref. 3, p. 107 ; Tennekes and Liunley, ref. 4, ch. 2) that this hypothesis 
does not have a firm physical base. Its most successful form, the mixing-length 
theory, requires a supplement of empirical information and theoretical modification 
(ref. 2, lecture 2), which diminishes the confidence with which it can be transferred 
from one type of flow to another. 

The question that is asked and explored in this report is whether an alternative 
zero- equation model of turbulence can be constructed by formulating a h 3 q)othesis for 
stress that is compatible with its known physical and mathematical properties and that 
retains simplicity of form by relying on a minimum of empirical constants and func- 
tions. Elementary observations are shown to lead to a simple expression for stress 
that approximates closely the distribution found in several common types of flow. The 
stress formula is applied to channel and pipe flows to find the velocity profiles and the 
variation of centerline velocity with pressure gradient. By using an approximate ve- 
locity profile in a system of integral equations, the friction coefficient, displacement 
thickness, and momentum-loss thickness in a flat-plate boundary layer are found. 
Comparisons are then made with experiment. 


HYPOTHESIS FOR TURBULENT SHEAR STRESS 

A hjrpothesis for turbulent stress should meet a number of requirements if it is to 
be a valid candidate for consideration and use in actual calculations. It should, first 
of all, lend itself to expression mathematically as a tensor inasmuch as the stress it- 
self is a tensor. If only one component of the stress (e. g. , the shear component) is 
being considered, it should be formulated in a way that could conceivably be generalized 
to a tensor of order 2. Second, it should be expressed as a function only of the dy- 
namic variable of the problem it is to be used in, namely, the time-mean velocity, its 
derivatives, and, possibly, its integrals. Finally, it should express the experimental 
fact (ref. 4) that the stress has a nonlocal character and that its value at any point de- 
pends not on the dynamic variable or its derivatives at that point alone, but on condi- 
tions at other points in the flow as well. 

As a restriction to problems of a simple yet practical sort, we confine our remarks 
to nearly unidirectional flows, steady and incompressible, whose lateral extent is much 
less than the extent in the streamwise direction. Among these are flows in pipes and 
channels, boundary layers, wakes, and jets. The velocity component in the direction 
of flow is U(Y). (Symbols are defined in appendix A. ) 

It is convenient to begin by considering the integrated frictional drag force on a 
body. A very elementary property of turbulent drag is that it is closely proportional 
to the square of the imdisturbed velocity of the stream. While it is true that in 
boundary- layer and pipe flows, for example, this behavior is modified by a weak 
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Reynolds number dependence, such deviation can be ascribed to the presence of viscosity 
rather than to the turbulent motions themselves. Not only is the square of the velocity 
descriptive of the drag in a quantitative way, it is also dimensionally appropriate when 
combined with the other available dimensional quantities, the density of the fluid and 
the characteristic dimension of the body. This suggests that at any point in the flow a 
component of stress may be related to the square of the local mean velocity. Further, 
if the shear stress, which is the component of interest here, were approximated by a 
polynomial of arbitrary power in the velocity, it would need to have no more than two 
roots since it vanishes only twice: once at the maximum mean velocity, and once at or 
near the minimum velocity in the flow. 

All of these considerations make it plausible to model the turbulent shear stress by 
a quadratic form in the local mean velocity, 

r = p(AU^ + BU + C) (1) 

If the two roots are located at the maximum velocity and at some ’’cutoff" velocity 
Ug below which the shear stress vanishes, then 

T = ap(U^ - U)(U - Uj,) (2) 

where a is a pure number that should be independent of mean-flow quantities. The 
quadratic form (eq. (2)) is, perhaps, the simplest h 3 ipothesis that could be made for the 
turbulent shear stress and has the additional property that it resembles the definition of 
Reynolds stress as a time-averaged product of fluctuating velocities by exhibiting a 
product of mean-velocity differences. It also has the desired appearance of a component 
of a second-order tensor, especially if oi were the component of a fourth-order tensor. 
In a very limited sense it is even a nonlocal description of the stress, introducing the 
velocities and at points other than that at which the stress is beii^ calculated. 

Actually, of course, any point in the flow is connected with almost all other points 
across the local cross section by the turbulent motions, so equation (2) should be thovight 
of as a first approximation to the most general description of this sort. 

The conic section described by equation (2) is a parabola in the mean velocity. 
Whether or not this is descriptive of the experimentally observed shear stress can be 
judged by comparing the hot-wire measurements of shear stress in flows of the type 
being considered here with parabolas fitted to them on the basis of least square error. 
This has been done in figure 1. The circular symbols represent actual hot-wire meas- 
urements or points from curves faired through hot-wire meauvu*ements. The sqviare 
symbols are values derived from mean-velocity data, either by estimating the cutoff 
velocity or by measuring the profile slope. The parabolas are required in each case to 
pass through the point of symmetry, U/Uj^ = 1, t = 0. 
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The agreement between the experimental data and the least- squares parabolas is 
sufficiently good that the error in using the latter would not be a large fractional part 
of the true value for flows boimded by walls, namely, bovmdary layers and pipe and 
channel flows. However, for free turbulent flows, the half jet and the wake, the agree- 
ment is not as good. The parabola is successful in representing the stress as obtained 
indirectly from mean-velocity measurements, but these data are not in agreement with 
the stress as measured directly by hot wire. The right leg of the parabola for the wake, 
which corresponds to the center of the wake, may yield a shear stress that is too small. 

This latter point is emphasized by wake measurements of Chevray and Kovasznay 
(ref. 5), which are shown in figure 2. (Data are presented in ref. 6. ) At the center of 
the wake, the slope of the data is nearly infinite and the parabola, which is chosen to 
pass through the points t = 0, (U^ - U)/(Ug - Uq) = 0, 1, cannot accurately represent 
the stress in the center region. (U^ and Uq are the velocities at the edge and center 
of the wake, respectively. ) 

On the basis of these observations the stress hypothesis (eq. (2)) will be used sub- 
sequently to solve the full boundary-value problem for only those flows that are bounded 
by walls, although in the following section it will be used to find similarity laws for free 
turbulent flows as well. 

One result of the comparison of hypothesis with experiment in figure 1 is the ap- 
pearance of a small but finite region of the flow, corresponding to 0 :< U :< U^., in which 
the shear stress vanishes. Although the hsqjothesis does not require that this thin layer 
of flow be steady and, therefore, laminar, it does state that in it, for all practical pur- 
poses, the random fluctuations of velocity are vmcorrelated. (However, the correlation 
of velocity and temperature fluctuations apparently does not always vanish over a finite 
region. This was shown by the calculations of Deissler (ref. 7) for high Prandtl num- 
ber. ) The appearance of a region of zero stress in this hypothesis corresponds to the 
influence of viscous damping functions in the eddy diffusivity of Deissler (ref. 7) and 
the mixing length of van Driest (ref. 8). 

The two adjustable parameters of the hypothesis are a and U^. The coefficient a 
occupies a position analogous to that of viscosity in viscous stress and should depend 
only upon the character of the turbulent motions; that is, it may be depend upon the 
geometry but not upon any dynamic parameter (imposed pressure gradient, free-stream 
velocity) that might be represented in the product of velocity differences. It may show 
dependence on preturbulence, free-stream turbulence, and so forth. 

The cutoff velocity, for flows boimded by walls, is in the region dominated by the 
skin friction and, in general, should be a function of the friction coefficient. The fimc- 
tion itself should change little from one bounded flow to another. In the calculations to 
be imdertaken subsequently, the functional dependence wiU be taken as the simplest 
kind. In figure 3 is shown a sketch of the velocity distribution near the wall in the fa- 
miliar u'*', y'*' coordinates. The cutoff velocity occurs at the jimcture of the profile of 
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the turbulent region with the line u'*' = y’*', which characterizes the siiblayer, in which 
only viscous stresses occur. It is assumed here that, for all flows of one class, the 
cutoff velocity u't = U /U will be a constant. The constant will be allowed to vary 

C C T , 

from class to class, but all pipe flows will have the same value of u^, all channel flows 
their own characteristic value, and so forth. 


SIMILARITY LAWS 

An immediate test that can be applied to the hypothesis (eq. (2)) is a comparison of 
the similarity laws that it predicts with those encountered e3q)erimentally. In deriving 
the similarity laws one assumes that the viscous stresses may be neglected compared 
with the turbulent stresses. The similarity laws predicted by the shear-stress hypoth- 
esis for plane and circular wakes, plane and circular jets, and the half jet are derived 
in appendix B. The results are shown in the following table: 


Flow 

Width 

Centerline ve- 
locity or ve- 
locity de icit 

Plane jet 

X 

X-I/2 

Circular jet 

X 

x-‘ 

Plane wake 

xV2 

X-1/2 

Circular wake 

xVs 

3J-2/3 

Half jet 

X 

x“ 


These are exactly the laws derived from the mixing-length theory and confirmed by ex- 
periment (Schlicting, ref. 9). 

The partial differential equation that governs the half jet is also appropriate to the 
outer portion of the flat-plate boundary layer. Consequently, the outer portion of that 
boimdary layer should have Y/X as a similarity variable for the velocity profiles. 

This was confirmed experimentally by Townsend (ref. 10). 


CHANNEL FLOW 

As a first application and test of the hypothesis for turbulent stress, consider fully 
developed flow between infinite parallel walls of spacing 2D (fig. 4). There is no ac- 
celeration of the fluid, so the pressure gradient and retarding forces are in balance. 
The maximum velocity occurs at the centerline, so we may write for the fluid in 
the upper half of the chaimel 



3X 


(3a) 


and 


BY 


II 0'P(U 

3Y 


m 


U)(U 



for U„ < U < 

c m 


i!u,5P torO<U<U, 

3y2 


(3b) 


where the lateral coordinate Y is measured from the channel centerline. 

If the reference velocity is chosen to be the maximum velocity that the fluid would 
have in laminar flow, namely, 


u .5! 

^ 2jtx 3X 

and if the half width D is chosen as the standard of length, equation (3a) after one inte- 
gration and equation (3b) after two integrations become 


“ - ap(Ujj^ - u)(u - u^) = -2y for < u < u^^ (4a) 


and 


u = 1 - y^ for 0 < u < Ug 


(4b) 


The new parameter, p, is the dimensionless pressure gradient or, equivalently, the 
Reynolds number based on the reference velocity. 



V 2ixv 


15 

3X 


(5) 


Equation (4b) states that the velocity in the sublayer follows the parabolic distribu- 
tion that it would have in a laminar flow at the same pressure gradient (ref. 3, fig. 9. 1). 
Equation (4a) is a Riccati equation (Goldstein and Braun, ref. 11) that can be trans- 
formed^ to the Airy equation in the following manner: The velocity deficit w = Uj^^ - u 
is replaced by a variable W through 


^This transformation was pointed out to the author by Dr. M. Goldstein. 
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w = - 


ap W 


There foEows a second-order equation 


W” + apw-W* + 2apyW = 
c 


which, in turn, is transformed by W = e®^T(y) to 

T” + (2s + apWg)T* + + apw^s 

with 


= u„ - u„ 

CSC 


The first derivative can be eliminated by choosing 


apw. 


s = - 


so that 


T»t + 


2apy - 


fCLWr 


T ^ 


Mroducing a new independent variable x defined by 


X 


where 


f 

\ 4 / 


a^(2o!p)2/^ 


leads to the Airy equation 


dfr 

dx^ 


xT = 0 


0 

h 2apyjT = 0 

( 6 ) 


0 


(7) 


( 8 ) 


(9) 
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In terms of the dependent variable T, the velocity deficit is 



where 


( 10 ) 


m 



( 11 ) 


is the value of x at the origin of y. 

The most general solution of the differential equation (9) is 

T = CjAi(x) + CgBiCx) (12) 

where Ai(x) and Bi(x) are Airy functions (Abramowitz and Stegun, ref. 12) of the first 
and second kinds, respectively, and and C 2 are constants. Because of the form of 
the expression (10) for the velocity deficit, only the ratio of the constants k s Cg/c^ 
enters into the solution: 


w = 2 Ter . Ai»(x) + kBiHx)' 
a Ai(x) + kBi(x) _ 


(13) 


The coordinate x has its origin at a point near the wall where the total stress 

(2y = TxD/ juU in dimensionless form) is equal to the maximum ap(w /2)^ that the tur- 

LX O 

bulent stress alone reaches in the flow. The coordinate of the cutoff velocity is obtained 
by combining equations (4b), (7), and (11), which yields 


At the center of the channel the velocity deficit vanishes to give 

^ _ Ar(x^) -f y^Ai(x J 


( 15 ) 


The boundary condition at x^ is 
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( 16 ) 




2 

a 



Ai»(Xg) + kBi'(Xj,) 
Ai(Xj,) + kBi(Xg) 


which, combined with equation (11) to eliminate w , yields a second expression for k, 


Ai^(Xc) -V^Ai(x^) 
Bi'(Xc) -V^Bi(x^) 


(17) 


Equating the two expressions for k relates the coordinate values x^ and x^ of the 
two velocities u^ and u^. The relation may be expressed in a simpler form by noting 
that estimates from experiment show 


= 0(10) 
\ = 0 ( 1 ) 


(18) 


Then, since for x > 2, Ai’(x), Ai(x) « Bi’(x), Bi(x), it follows from the first expres- 
sion for k (eq. (15)) that 


k « 1 (19) 

and from the second expression for k (eq. (17)) that the relation between x^^ and x^ 
is approximately 



Ai’(x^.) 

Ai(x^) 


( 20 ) 


Introducing equations (11) for x^ and (14) for x^, supplemented by equation (8) for a, 
yields a relation among u^^, u^, and p. The cutoff velocity can be eliminated because 
we have set u^ = Constant and, by definition. 



c r 


= u. 


Ur Ur 



or 
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( 21 ) 


u 


c 



The resulting dependence of centerline velocity on the pressure gradient p is 
plotted in figure 5, where it is compared with the ejqjerimental data of Laufer (ref. 13). 
The values of a and u^ that were found by trial to provide a good fit were 0. 01 and 6, 
respectively. The agreement between theory and experiment is good, the maximum dif- 
ference being at the high values of pressure gradient, where it is about 6 percent. 

The dimensionless average velocity through the channel, defined as 


r 


U dY 


UrD 



can be found from equations (13) and (4b) to be 

- + "r 2 -Ai(Xf.) + kBi(x ) 2 - y (2 + u ) 

u =-i2 5 y +_^ln ^ ^ 2 ^ 

2 "" 3/2 Ai(x )+kBi(x ) 3 


where 


- "c 

This quantity is also shown in figure 5. The ratio approaches 1 as the pressure 

gradient becomes very large. 

In figure 6 a velocity profile at a Re 3 molds number (Re = pu^^) corresponding to one 
of those at which Laufer *s measurements were taken is compared with the experimental 
profile. The theoretical profile is not as full as the experimental profile, and this must 
be ascribed to a failure of the stress hypothesis to provide a sufficiently strong stress 
in the center of the channel. Nevertheless, because the centerline velocity is predicted 
accurately, the mean velocity is not likely to be in substantial error. 


PIPE FLOW 


The momentum equation for fully developed pipe flow is (Laufer, ref. 14) 

1 dP _ 1 d ^ ^./d^U ^ 1 du\ 
p dX R dR p \(Jr 2 R dR/ 


( 22 ) 
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where, as usual, X and R are the axial and radial coordinates, respectively, Ihtro- 
ducii^ equation (2) and nondimensionalizing the radius of the pipe Rq and the centerline 
velocity for laminar flow 


Kq 

Ur= — 
^ 4/n 


dX 


(23) 


yield 




where 


PQ " 


V 4fJLU 


dP 

dX 


r 



(24) 


(25) 

(26) 


The reference velocity (eq. (23)) is just one-half that for a channel flow with its 
half -width equal to the pipe radius; the factor 1/2 is also reflected in the definition of 
the dimensionless pressure gradient (eq. (25)). 

Equation (24) governing the flow is identical to the Riccati equation treated for chan- 
nel flow, and the solution proceeds in exactly the same way. Figure 7 shows the de- 
pendence of the dimensionless centerline velocity upon the dimensionless pressure 
gradient Pq as compared with the experimental data of Laufer (ref. 14) and Sandborn 
(ref. 15). The value of the parameter a has been kept at ^.01, but the value of u^ 
has been raised from 6 to 7 to give the best fit to the experimental data. The calculated 
curve is slightly above the data at low p and slightly below at high p. Further calcula- 
tions show that raising u^ above 7 can place the curve through the high-p data and that 
lowering it can place the curve through the low-p data. It can be concluded, therefore, 
that the choice of constant n't is nearly the correct one but that improvement could be 
obtained with a more general dependence on the parameter p. 

Also in figure 7 the dimensionless average velocity 

u = 2 ^ ru dr 
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is compared with the universal resistance law for pipes (Schlicting, ref. 9). The cal- 
culated u line is nearly parallel to the u^ line, which shows that the calculated ve- 
locity profiles do not change with p sufficiently for the mean velocity to follow the ob- 
served values closely. 

Figure 8(a) shows the velocity profiles at Re = of 250 000 and 25 000. Just 

as in the case of channel flow, the calculated profile has the correct overall character- 
istics but is not as full as the profile found experimentally by Laufer. The same pro- 
files are also given in terms of the similarity variables u"^ and y"*" in figure 8(b). This 
plot shows both of the effects already observed in figures 7 and 8(a). The calculated pro- 
files tend not to be as full as the observed data, and the centerline velocity is too high 
at the low dimensionless pressure gradient and too low at high p. 


FLAT- PLATE BOUNDARY LAYER 
Differential Equations 

For a boundary layer the hypothesis for the turbulent stress has the form 

T = ap(U 5 - U)(U - U^) for < U < Ug (27) 


where Ug is the stream velocity at the edge of the boundary layer. We again assume 
that the minimum velocity in the turbulent region is given by 



Uc 

— = Constant 
U_ 


(28) 


with customary way. 

The stress hypothesis (eq. (27)) is introduced into PrandtPs equation for the incom- 
pressible boundary layer 


U 


ax 


3Y 




3 T 
3Yp 


(29) 


and at the same time the dependent and independent variables undergo the transformation 
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( 30 ) 


U = Ugu, V = aUgV 
X = Lx, Y = Q!Ly 

where L is a typical length for the geometry and u, v, x, and y are considered to 
be 0(1). 

Since a = 0. 01 approximately, judging from experience with the pipe and chaimel, 
transformation (3) states that a boundary layer should be about 100 times as long as it is 
thick. This is in accordance with observation. Prandtl’s equation is now 


3u 0u 1 ^^6 /, 2\ 1 

u — + V — = (1 - u ) -i 


9x 


0y Ug dx 


— (1 - u)(u - u^) 


Re 0y^ 


for u. < u < 1 
c 


(31) 


where a local Rejmolds number 


Re 

V 


(32) 


has been introduced. In this form of the equation, the inertia terms and the turbulent 
stress terms are of the same magnitude, but the viscous term has a small coefficient. 
This is in keeping with the well-known property of turbulent boundary layers that the 
viscous forces are appreciable near the surface, where the velocity derivative is high, 
although the turbulent stresses dominate throughout most of the boundary layer. 

The continuity equation 0U/0X + 0V/0Y = 0 becomes under the transformation 
(eq. (30)) 


9u ^ 0v 
9x 0y 


u 



dx 



(33) 


It will be convenient subsequently to have the properties of the velocity profile near the 
wall in terms of a and Re. Since the wall stress is given by = (jitdU/dY)Q, the 
dimentionless velocity derivative at the wall can be expressed in terms of the friction 
coefficient f s 2 t^/ pUg by 



= 0 ! Re — 
2 
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Moreover, Pohlhausen's compatibility condition at the wall (Watz, ref. 16, p, 80), found 
by evaluating the motion with u = v = 0, requires 


1 ^^5 , 1 

ye 3y2 


Therefore, near the wall, 


u 


_ O' Re f y _ d Re y2 ^ 
2 2 dx 


(34) 


Moments of the Equation of Motion 

Solutions to the partial differential equations governing the boundary layer have been 
found principally in three ways. For a laminar boundary layer, it has been possible to 
reduce the system to one ordinary differential equation by a similarity analysis, at 
least for certain pressure distributions. This is not a viable approach for the system 
described by equations (31) and (33) because in the turbulent boundary layer the viscous 
stress and the turbulent stress lead to different similarity laws and it is not possible to 
arrive at an ordinary differential equation. A very direct method, which is as available 
to the turbulent boundary layer as to the laminar, is numerical integration by a march- 
ing technique. It would be necessary to use a smaller grid size in the viscous layer 
than in the outer, turbulent layer and, in addition, it would be necessary to follow ac- 
curately the boundary between the two regions at u = u . Thus, the marching tech- 
nique as now practiced for laminar boundary layers, although requiring some modifica- 
tion, is adaptable to the turbulent boundary layer as formulated herein. 

The method of solution that will be used here is the integral moment method, in 
which the equation of motion is multiplied by a weighting function and then integrated 
over the boundary-layer width to produce an ordinary differential equation in the inde- 
pendent variable x. Accounts of its development and practice are given by Walz 
(ref. 16) and Tetervin and Lin (ref. 17). Most frequently, the first velocity moment is 
used to provide an energy equation featuring a dissipation integral. However, Tetervin 
and Lin derived the most general form of an integral moment equation by usii^ arbi- 
trary powers of u and y, and this generality makes it possible to choose a preferred 
set of the equations. 

The present, brief experience with the pair of equations (31) and (33) has led to a 
preference for a system of moments in (6 - y)”, that is, taking moments about the edge 
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of the boundary layer. A reason for choosing coordinate moments over velocity mo- 
ments may be seen by considering the zeroth moment of equation (31), which is (ref. 16, 
p. 87) 


+ _L ^ (6* + 0) = X. • (35) 

dx Ug dx 2a. 

where 9 and 6* are the dimensionless momentum-loss and displacement thicknesses. 
This equation, the lowest order equation of any moment system, has three dependent 
variables, 0, 6*, and f. If the next moment equation to be added is the first velocity 
moment, a new d 3 mamic variable, the energy thickness, is introduced as an imknown. 
Thus, from this point of view the system does not close: the number of unknowns al- 
ways exceeds the number of equations. It seems more appropriate to first add two 
y-moment equations to equation (35) so that the three equations may be thought of as a 
system for 0, 6*, and f and only then, if another equation is required, to add the first 
velocity-moment equation and the energy thickness to the system. 

Moreover, when a system of velocity-moment equations was used to solve the pair 
of equations (31) and (33), it was found that the system was \mstable. This may have 
been due to the choice of approximate velocity profiles, which have to accompany the 
moment system, but it is more likely that the cause lies in the high powers of u, which 
multiply the error for any assumed profile and cause it to grow rather than diminish. 

The reason for choosing 6 - y rather than y as the moment arin is one of con- 
venience. A moment in (6 - y)^ is no more than a linear combination of moments in 
y*^, y, . . .,y”~^, y^ and, hence, introduces no new information. However, since mo- 
ments in y tend to yield more complicated equations than moments in u, any simpli- 
fication that can be attained is desirable. The derivation of the moment equations in 
6 - y is given in appendix C. It is apparent there that a number of terms arising from 
partial integration vanish at the outer edge of the boundary layer due to the choice of 
6 - y for the moment arm. In addition, as a subsequent discussion points out, some 
suggestions for velocity profiles are expressed in terms of 6 - y and, as a result, in- 
tegrals of moments in 6 - y are simplified. 

As shown by equation (C6) of appendix C, the n^^ integral moment equation for mo- 
ments taken around the edge of the boundary layer and for the restricted case of the flat- 
plate boundary layer is 
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uv dy 


t V 0 f <ly f + V „,(«) * n (5 - y)»- ' 


6^— ) + 6 . + n(n - 1) /" (6 - y)" ^ u dy 

Vay/Q ’ -^0 


“ f (® • y)”"^ (1 - u)(u - u )dy 


( 36 ) 


a" Re 


where 6 . is the Kronecker delta. 

J 

We shall require explicitly only the first three of these equations, for n = 0, 1, 
and 2. These are 


For n = 0: 


For n = 1: 



d s 2 , d6 

— / (6 - y)u dy - — 

dx ‘>'0 dx 




uv dy 


(37a) 


1-6 


9u 

9y. 


Re 


rb 

J (1 - u)(u - u^)dy 


(37b) 


For n = 2: 


— f (S - y)^u^ dy - 2 — /* (6 - y)u^ dy + 2 /* (6 - y)uv dy 
dx *^0 


dx *^0 


Re 


u dy - 6YiH 
‘'O 


3yy 


rb 

2 y (5 - y)(l - u)(u - Ug)dy 


(37c) 
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The zeroth equation is just a special form of equation (35) for a constant stream 
velocity. It is obtained after writing, in equation (36), 




.a d6 
u dy + — 

dx 


The zeroth equation is noteworthy for being independent of the form assumed for the 
turbulent stress; only the viscous stress at the wall enters into it. The n = 1 equation, 
however, contains an integral of the turbulent stress across the tmrbulent portion of the 
boundary layer. This integral represents the total force exerted by the plate laterally 
against the fluid. The equation for n = 2 has an integral of the first moment of the 
turbulent stress; and since the stress peaks close to the wall, this shows that taking the 
moment about the edge of the boundary layer is once more desirable. It emphasizes the 
effect of the stress more than if the moment were taken around y = 0. 


Choice of Velocity Profiles 


The integral moment method proceeds by performing the integrals and derivatives 
in the system (37) on a velocity profile that contains as many free parameters as there 
are equations in the system. The profile must have the form of equation (34) near the 
wall; and it is desirable that at the edge, y = 6, of the boundary layer, 3 u/8y = 0 in 
agreement with equation (31). 

Power-law profile. - One of the most successful representations of the dimension- 
less velocity has been the form u=(y/6)^/®, where e > 1. Unfortimately, it does not 
meet the requirements set previously for behavior at the wall and the edge of the bound- 
ary layer. The failure at the latter point is not serious; but in order to avoid a singular 
slope at the wall, we use equation (34) (modified for constant Re) near the wall so that 


u = g..%. ^ y 
2 


for 0 < u < Ug 


(38a) 


u = 



for u^ < u < 1 


(38b) 


If we again choose u^ by requiring that the bomidary between the sublayer and the 
timbulent region be at a constant velocity on the u"*", y"*" diagram, we find that 
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U U;: 
c 0 





or 



It proves convenient hereinafter to use yjt rather than f or to choose as a new vari- 
able 


F s 



(39) 


Then 


^c 


J 

F 


(40) 


where 


Js 



= Constant 


a Re 


(41) 


The profile (eq. (38a)) becomes 


u 



(42) 


The two parts of the profile (eqs. (38)) must be matched at u = u^, leading to a 
connection between the three free parameters 6, e, and F of the profile. 


6 = 


J 


pl+e 


(43a) 


Its differential form is 
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(43b) 


— + (In F)e + (e + 1) = 0 

6 F 


where 6 s d6/dx, and so forth. 

In addition to using this relation between the parameters, we shall also use the 
first two of the moment equations (37). In appendix D it is shown that they yield a sys- 
tem 


^01^ + ^02® + ^03^ = ®0 

(44a) 

Aii6 + Ai26 + A^3F = Bi 

(44b) 


and we bring equation (43b) into this system by writii^ it as 

^ 31 ® + ^ 32 ® + ^ 33 ^ = ®3 

with A.^-^ = 1/6, Ag 2 = hi F, Agg = (e + 1)/F, and Bg = 0. 

Initial conditions . - The initial conditions for the system (44), however originally 
specified, must be converted to the parameters 6, e, and F. To illustrate the integral 
moment methods we shall integrate the flat-plate, boundary-layer system (44) by starting 
at station 4 of the Wieghardt experiment and using the measurements as given in the re- 
port of the Stanford Conference (ref. 18) for starting conditions. At X = 0. 387 meter, 
the dimensional momentum-loss thickness © = oLd is 0. 0924 centimeter and the friction 
coefficient f is 0.00364, where L s 5 meters is the length of the plate. 

The initial condition for F comes from equation (39). According to appendix D, 




+ 6e^ 


1+ F' 


e+1 


e + 1 


1 - F^^ \ 
e+ 2 / 


(45a) 


which, combined with equation (43a), yields 

1 F ^ e /l - F^^ 1 - F^^ \ 0 

2 3 F®+1\ e+1 e + 2 j aLJ 


(45b) 


as an initial condition for e. The remaining initial condition for 6 follows from equa- 
tion (43a). 
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Three-parameter profile. - A velocity profile introduced by Pai (refs. 19 and 20) 
for channel and pipe flows and later adapted by Sandborn (ref. 21) to boundary layers 
has the form 


^ = (46, 

This is a profile with five parameters that is intended to hold over the entire boundary 
layer. As in the previous example we shall use this form only in the turbulent region 
and supplement it with a straight-line profile in the sublayer. We shall also change the 
square term to a cubic in order that the second derivative vanish at the edge of the 
boundary layer as required by the- equation of motion. To meet the boundary condition 
at the outer edge of the layer, we set A equal to 1; and to reduce the number of pa- 
rameters, we match both the values and slopes of the two profile segments at u = u^. 

The velocity profile is then specified by 

u 

u = — ^ y for 0 < u < u (47) 

yc 

where u = F and y = J/F as before; and 
c 

u = 1 - b?^ - c?^® for Up < u < 1 (48) 


where 




Lzjl 


6-y 


c 


(5 - yjF" 


c = 


3J 


- 1 + F 


e - 1 


b = 1 - F - c 


Since the matching conditions have determined b and c, only 6, e, and F remain to 
be determined by the moment equations. Using the moments specified by n of 0, 1, 
and 2, we again get a system of ordinary differential equations of the form (44) whose 
coefficients are defined in appendix E. 
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Calculations 


The two systems of moment equations corresponding to the two approximate veloc- 
ity profiles just described have been integrated, and the results are shown in figure 9. 
The two-parameter profile can represent the momentum-loss thickness well but, since 
the two initial conditions were used for 9 and f, the displacement-thickness curve 
shows an error of about 8 percent at the end of the plate. It is apparent that the solu- 
tions want to give the flow a constant friction coefficient. This can probably be attri- 
buted to the very simple assumption made in choosing the cutoff velocity. A weak de- 
pendence of u^ upon f or F might improve the agreement with experiment at the ex- 
pense, of course, of adding an additional experimental parameter. 

The three-parameter velocity profile shows some improvement over the two- 
parameter profile in the representation of displacement thickness but is hardly otherwise 
superior. 

Figure 10 compares the two velocity profiles with the measured velocity distribu- 
tion. The velocity measurements and the calculated profiles are plotted against the 
similarity variable for the outer portion of the boundary layer = Y/(X - Xq), where 
the virtual origin has been fotind by trial to be Xq = -27 centimeters. Both profiles 
have difficulty in duplicating the knee of the curve in the region near the wall. In the 
case of the power-law profile this can be attributed to the power itself, which instead of 
ranging from 7 to 9 is approximately 4. Presumably, this is due to the introduction of 
the straight-line segment in the sublayer. The three-parameter profile follows the cubic 
term from y = 6 to nearly y = y_ because the value of the exponent e is so large. At 
that point the third term provides a sharp swing from the cubic curve down to the point 
u = Up, 7} = at the end of the linear profile of the sublayer. For each case the point 
Uc, Tjp is marked by a heavy dot. 

The adjustable parameters a and u^ turn out to be somewhat profile dependent in 
the case of the boundary layer. For the two- parameter profile the values a = 0. 01, 
u]^ = 6, which were used for the channel, were again employed. However, the values 
that prove appropriate for the three-parameter profile are a = 0.0135 and u^ = 10. 


CONCLUDING REMARKS 

Comparii^ the calculations with experiment indicates that the stress hypothesis has 
the proper overall characteristics for describing turbulent flows near walls. The gen- 
erally correct trends that the theory exhibits are probably attributable to the appearance 
of the square of the mean velocity in the stress; the lack of precision over a wide range 
of operating parameters (p or x) lies in the choice of the cutoff velocity. It appears 
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that the dependence of the cutoff velocity upon the friction coefficient should be refined. 

Nevertheless, the variation of cutoff velocity with friction coefficient in the bound- 
ary layer, although not precise, is encouraging for the description of the flow near 
transition and under accelerated free streams. Because the cutoff velocity increases as 
the square root of the friction coefficient, the theoretical boxmdary layer will tend to 
become mostly sublayer whenever the friction is high. This is just the condition ob- 
served ejq)erimentally near transition and relaminarization. An accurate description 
under such conditions will require a more flexible profile in the sublayer. 

It is not apparent whether the argument that led to the proposed expression for the 
shear component of stress can be extended to provide the normal stresses as well. Al- 
though the normal stresses do not arise in the particular examples discussed herein, it 
is agreed that they may be important in some flows, in a boundary layer near separation 
for example, and certainly will be needed for the calculation of pressure distributions. 

In addition, a representation of other components of the stress may give some clue to 
the failure of the h 3 q)othesis in free turbulent flows. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 20, 1977, 

506-24. 
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APPENDIX A 
SYMBOLS 


A,B, C 

Ai(x), Bi(x) 
Ai*, Bi’ 
a 


b, c 


D 


e 

F 

f 

J 

k 

L 


n 

P 

P 

PO 

R 

Rq 

Re 

r 

Si 

s 


constants 

matrix coefficients 
Airy fimctions 

derivatives of Airy fxmctions 

(2c*p)2/3 

component of column vector 

coefficients in velocity profile 

function of 6, e, and F defined in appendix E 

constant 

half-width of channel 
form parameter 

friction coefficient 

(<) 7 “ 

constant, Cg/cj 
length of flat plate 
order of moment 
pressure 

dimensionless pressure gradient, U^D/u 
dimensionless pressure gradient, 
radial coordinate 
radius of pipe 
Reynolds number 

dimensionless radial coordinate, R/Rq 
function of 6, e, and F defined in appendix E 
parameter in channel solution 
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T 

U 

u„ 


u. 


m 


U£ 

u. 

Ur 

U 


R 


w 

c 

u 

V 

V 

W 


w 

X,Y 

Xn 


x,y 


X 



a 


A* 

6 


6 * 


0 
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intermediate variable in channel solution 
mean velocity in X-direction 
cutoff velocity 
maximmn velocity in flow 
velocity at edge of boundary layer 
reference velocity for pipe, Rq | 3P/ 3X | /4/li 
reference velocity for channel, j 3P/ 3X | /2p 

friction velocity, 
dimensionless velocities 

C' T 

average dimensionless velocity 
velocity in Y-direction 
dimensionless velocity in y-direction 
intermediate variable 
dimensionless velocity deficit, u^ - u 
dimensionless velocity deficit, u^^ - u^ 

Cartesian coordinates 

virtual origin of boundary layer 

dimensionless Cartesian coordinates 

independent variable for Airy equation 

value of X at Y = 

value of X at y = 0 

value of y at u = u^ 

variables defined in appendix D 

adjustable parameter in stress hypothesis 

displacement thickness 

dimensionless nominal thickness of boundary layer 
dimensionless displacement thickness 
(6 - y)/(6 - y^) 
momentiim-loss thickness 


■ 


d dimensionless momentum-loss thickness 

ju. viscosity 

V kinematic viscosity 

p density 

T turbulent shear stress 

sum of viscous and turbulent shear stresses 
shear stress at wall 
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APPENDIX B 


SIMILARITY LAWS OF JETS AND WAKES. 
The equations governing jets and wakes are 



— (Y^U) + — (yV) = 0 (B2) 

ax BY 


where j = 0 for plane flows and j = 1 for axisymmetric flows. 


Jets 

Integrating the momentum equation (Bl) over the appropriate interval after first 
multiplying by Y^ yields 



Y% dY + 

ax 



yjy dY = 
aY 



where b = -°° for j = 0 and b = 0 for j = 1. The second term on the left can be in- 
tegrated by parts using equation (B2) to obtain 



yV — dY = 
aY 




Y^U— dY 

ax 


Since t and V vanish at 


the limits of integration, 


the momentmn integral equals zero. 


_a_ 

ax 



Y^U^ dY = 0 


(B3) 


We now introduce a similarity transformation 
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u = UpX f(?7), X = — y = r] = ^ (B4) 

where Uj^ is a reference velocity. We may take the first term in equation (Bl) as t 3 npi- 
cal of inertia terms and find it to be 


uiU = % x^“'^f(mf - nrjf) 

ax V 


The stress term becomes 


(B5) 


yj 9Y I p) 


c>u| ^2m-n ^ 
V ^3 dr] 


■ '>(* - y] 


(B6) 


where f and f„ correspond to and U and we have assumed that f„ is inde- 
m c ^ m c c 

pendent of x. 

Requiring equal powers of x in the inertia and stress terms (eqs. (B5) and (B6), 
respectively) yields 


2m - 1 = 2m - n 


(B7) 


or n = 1. Thus, the similarity variable for both plane and circular jets is t] = y/x. 
The decay law for the centerline velocity is given by equation (B3); 


i-x2m+(j+l)n 
3x “'b 

The power of x in this expression must vanish under the constraint (B7) so that 

m = - 

2 

or m = -1/2 for plane jets and m = -1 for circular jets. 


Wakes 

In this case we introduce the velocity deficit U’ = Um - U and consider only the 
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wake far downstream, where U* is small compared with The similarity trans- 

formation (eq. (B4)) is again employed, with U replaced by U’. The typical inertia 
term becomes 


2 

u — = - x“^"^(mf - n? 7 f*) 

ax V 


and the stress term becomes 


Q! ^ [T]jf(f - f)] 

where it has been necessary to assume again that can be considered constant. 
Equating powers of x in the inertia and stress terms yields m - n = -1. 
Condition (B3) becomes for the wake 

_a_^m+(j+l)n /'“°^jfdr, = 0 
0x *'b 


from which m + (j + l)n = 0. The solutions of the conditions for m and n are 


m = 




for plane wakes 


m 


p 

3 


> 



for circular wakes 


Half Jet or Mixing Layer 

The equations of motion for a half jet or mixing layer are those for a plane jet. We 
take the lower fluid to be at rest and the upper fluid to have the constant, undisturbed 


28 



velocity The latter requirement forces m to equal 0 in equation (B4). Intro- 

ducing this result into the comparison of inertia and stress terms in the plane- jet 
analysis yields n = 1. 



APPENDIX C 


moment equation for flat- plate boundary layer 

The integral of the n moment by 6 - y of the dimensionless boundary- layer equa- 
tion (31) with Ug constant is 


\ dx ay/ 


(6 - y)' 



1 au 

dy 

Re 


It is desired to transform this equation to a manifestly ordinary differential equation in 
X. Beginning with the second term on the left side, one has 


/ (6-y)”viSdy= /■ (6-y)”iMdy 

^ 9y ^ ay 


-/ 

^0 


(e-rt-uiidy 


J [(6 - y)” uv]dy + n J uv(6 - y)”"^ 


+ / (6-y)''uiHdy 


= nv(6) + n /* (6 - y)”"^ uv dy + /* (6 - y)“ u — dy (C2) 


The last term in equation (C2) is just equal to the first term in equation (Cl) and their 
sum is 


(6-y)"iH-dy=-l / (6 - y)“ dy . 

dx dx / 


6 /*^ 

[(6 - y)” u2]^ 6 - n6 / 


(6-y)“"^u^dy 


— A (6 - y)“ u^ dy - 5 0® ■ <*y (C3) 
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On the right side of equation (Cl) the first terna is proportional to 

/’°(6-y)”i!sd,= r(6-y)"isl%„ / iS (6 - y)”'! dy 

= - + 6 j + n(n - 1) f u(6 - y)""^ dy 

Wo ^ -^0 


and the second term is 


- y)”^ (1 - - V*^y = - y)”(i - - ^c>]\ 


n y' (6 - y)”"^(l - u)(u - u^,)dy 


= n y* (6 - y)"‘^(l - u)(u - u^.)dy 


Combining equations (C2) to (C5) yields for the entire equation 


(5 - y)^ dy - 6^^ q6 - n6 (6 - y)”"^ u^ dy+ 6^^^ qV(6)+ n (6 - y)”^ ^ uv dy 


[ ^4 I 


(6 - y)""2 u dy ^ 


2 „ 
a Re 


+ n /* (6 - W (1 - u)(u - u^)dy (C6) 

^ O 


Now consider, further, the case when the velocity profile is specified as a functiofn 
of three parameters, one of which is the nominal boundary- layer thickness 6: 


I 



u = u(y; 6, e, F) (C7) 

We can transform equation (C6) to an ordinary differential equation in 6, e, and F by 
introducing from the continuity equation 


v(y) = - — dy = - (— 6 + — e + — p'^dy 

^ 3x \36 3e 3F / 


and 


v(6) = - — C udy+6 = - fs — + e — + F — 1 // u dy + 6 

dx ^ L 36 3e 3FJ 

• • • • 

where 6 s d6/dx, and so forth. The coefficients of 6, e, and F on the left side of 

equation (C6) are then, respectively, 


\l = 


36 L 


/ 5 

(6 - y)^'^ dy - 6 


3 


dy 


-n /\6-y)"-lu r^dr 

Jq Jq 


dy (C8a) 




dy 

(C8b) 


A 3 = — / (6 - y)“ dy - 6 Q J- / u dy 

3F ^ ^ 3F ^ 


- n (6 - y)"'l u r ^ dy’ dy n = 0, 1, 2 

Jq Jq 


(C8c) 


With the further definition for the right side of equation (C6), namely. 
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®n = — 


Re 


®nl-5 


” f—\ + n(n - 1) f (6 - y)“"^ u dyl 

\3y/Q Jo J 


n (6 - y)“"^(l - u)(u - u^)dy (C8d) 


we can rewrite equation (C6) in the desired form of an explicitly ordinary differential 
equation 

^nl® + %2® + An3^ = Bn 


In principle, any set of three of these equations may be used to solve for 6, e, 
and F. In practice, the simplest to use will be the first three: n = 0, 1, 2. 
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APPENDIX D 


MATRIX COEFFICIENTS AND INITIAL CONDITIONS FOR 
POWER- LAW VELOCITY PROFILE 

We assume the velocity is described by 


u = — y for 0 < u < u 


( 42 ) 


u = 


1/e 


for u„ < u < 1 
c 


(Dl) 


The coordinate at the partitioning velocity is 


^ F 


( 40 ) 


and the relation between the parameters 6, e, and F resulting from the matching at 


U = U^ IS 


6 =. 


4+e 


( 43 ) 


The dimensionless displacement thickness 


A* 

aL 


= (l-u)dy = ^‘"^l-^y^dy+ 1- 


1/e 


dy 


e + 1 \ pl+e/ 2 


Similarly, the dimensionless momentum thickness is 

•6 


0 

aL 


H 0 = / u(l - u)dy = - 1 + -J_r_e_ (1 . F®+^) - (1 - F®+2)]1 

•'O 1^2 3 „l+e Le + 1 e + 2 Jj 
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The coefficients appearing in the first two integral moment equations are obtained by 
substituting the velocity specifications (42) and (Dl) into equation (C8), performing the 
integrals, replacing y with equation (40), and then performing the differentiations with 
respect to 6, e, and F with J held constant. In terms of the definitions z = y^/S = 
J/F5, Z 2 = and Zg = 1 - Zp the coefficients are 


^01 = ^1 ■ ^2 + 777 - 77 “ 

e + 1 2 + e 


^02 ~ ® 


1 - Zj^ Zj^ In z 2(1 - Z 2 ) 2 z 2 in z 


L(e+1)2 e(e+l) (2 + e)^ e(2 + e)J 


"7 ■ ^ 2 ^ ■ f 


^0 = 


2 t t , 

a J Re 


1 - Zo 1 “ z- 


= 6e«< 


e + 2 e + 1 


1 - z 


2(e + 1) 


^ - (1 H- Zj) 


i ^2 - 4 ) 


Ai2 " 5' 


2(1 - Zj) 1 - zj 


2 -1 


,(e + 2)^ 2(e + 1)^ 


+ 6 


21 2^2 


^1 \lnz 


1 2 + e 1 + e/ e 


(e + ir L 


z^ — i Jin z + — ^ (2 + Zo) — 
^ 2 / e+ 1 ■^4 




J e6(l - zj) 

4 e + 1 


-6 

Q!^ Re \ ^ / 


e(l - Z 2 ) e(l + F)(l - z^) 


2 + e 


e + 1 


+ F(1 - z) 


A third differential relation between 6, e, and F is obtained by differentiating 
equation (43). After multiplication by 6, the coefficients are 
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“^21 “ ^ 

A22 = 6 In F 

A23 = — (e + 1) 

B2 = 0 



APPENDIX E 


MATRIX COEFFICIENTS AND INITIAL CONDITIONS FOR THE 
THREE -PARAMETER PROFILE 

If the velocity profile defined by equations (47) and (48) is used to evaluate the 
coefficients defined in equation (C8), they can be written as follows: 


Aqi = Sj + (6 - y^) 


F^S, 


3J(e - 1) 


Aq 2 = c(6 - yj (3S 


c^ r“3 


e - ly 


J Js. 

^03 = - 3 + ^ - yc>(®2‘'F - S4) 

F 


®0 = 


1 F^ 
a^Re ^ 


^11 = - yc>[2C7 - Cg - Sj + (1 - F)C^ + C^ J + (6 - y^.)^C5Cc 


Ai3 - J 


6-v 


6 

3'*‘ J.2 


^12 = - yc^‘ 

^ (Cj^i^ + 2Ci^) 


3cCg l2 + 


1 \ 3c(C^ 

3e + 1/ /o^ , 


- CJ 


(3e + 1)^ 






F - y^(6 - y^)Ci + (6-y^)24c2-^i-^+CpC 


Bi = (l - t (6 - y^)[(l - DSji - S 12 ] 

Q!^ Re \ '^ / 


A 21 = (5 - ye)^[Cio - 1 + 2(1 - F)Ci 4 + 2C^i - 2Cig] 





2c 


(3e + 1 )^ L 


3C 


(3e + l)Cj^g 


14 


- C 


11 


■*■ ‘^e^l6| 


a^3 = 4^-'^4V <= - ^ <= - yc)^^-^-FCxe, 


,3 2 5 


'c' •'c 


12 2 


^2 = 


Re 


62f2 


+ J + 2(6 - y^)C^ 


+ 2(6 - y^)2[(l - f)C7i - C^g] 


In these expressions for the matrix coefficients the following groupings have been used: 
y^ = J/F, c = [(6 - yj.)F2/3J - 1 + F]/(e - 1), and b = 1 - F - c 


a _ b , C 

Sii - _ + 


4 3e + 1 


o _ b2 2bc . c2 
®12 


7 3e + 4 6e + 1 


®1 ■ ®11 " ®12 


®21 “ 


1 1 


3e+ 1 4 


^ \7 3e+4 6e+l/ 


83 = i 


(3e + 1)^ (3e + 4)2 (6e + 1)^ 


S =l-2(^ + —^) 

^4 \7 3e + 4/ 


= 1 - ^ 


4 3e+ 1 


P _ 1 b c 
^2 


5 8 3e + 5 
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^3 = 


(3e + 2)^ (3e + 5)^ (6e + 2)^ 


C4 = 


3© + 2 3© + 5 6© + 2 


c -i- 1 

^5 “ 


4 3e + 1 


^6 " ^ 5^1 ■'■ 


C4 Cg 


3e + 1 4 


P _ b c 

C^l _ _ + 


5 3e+ 2 


^ . 2bc . 

^.rj2 + + 


8 3e + 5 6e + 2 




Cg = 1 + 2 (C^ - 1 ) + S 


12 


Cg = 2C2 - 2C4 - Cg 


CjQ=l-b--^ + ^ + J^ + -i 


e+1 3 e+2 2e+l 


Cii - 


e+1 e+2 2e+l 


r. _ 1 2b 2c 
*^10 - ■*■ 

3 e+2 


^13 ■ 


(e + 1)2 (e + 2)2 (2e + 1)2 
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r. _ 1 b C 




2 5 3e + 2 


^15 = ^5^14 


1 ^11 ^12 


3e + 1 3 24 


_ Cj2 - 2Cii 

^16 ; ■^'^15 


'191 


'192 


ec 


3e + 1 


^19 " ^191 ^192 


- 3(-CjgCj + CjgjCg + Cj^92^4) 


^18 " ^^ 19^14 ■ 2 ^ 191^12 ‘ ^ 192^11 


‘^6 = 


3J(e - 1) 


^e=7T7 




e - 1 


If A* and © are the initial values of the dimensional displacement and momentum- 
loss thicknesses, and if f is the initial value of the friction coefficient, the correspond- 
ing initial values of the parameters 6, e, and F are obtained by calculating displace- 
ment thickness, momentum-loss thickness, and friction coefficient from the profiles 
in equation (19) and equating them to the experimental values. The result is a set of 
three equations 
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- yc^Si2 = 


A* - e 

OfL 


for 6, e, and F. 
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(a) Boundary layer, dP/dx = 0. 
(From ref. 22. ) 




(d) Boundary layer, dP/dx>0. 
(From ref. 24. ) 


(e) Pipe, Re= 25 000. (From 
ref. 14.) 



(b) Channel, Re = 30 800. (From 
ref. 13.) 


(c) Half jet. (From ref. 23.) 



(f) Plane wake, Re-680. (From 
ref. 3. ) 


Ratio of mean velocity to maximum velocity, 

Figure 1. - Turbulent shear-stress distributions (varying scales). 
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Figure 2. - Shear-stress distributions in wake of a plate. 
(From refs. 5 and 6. ) 
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Figure 6. - Velocity distribution across a channel, a = 0. 01; uj * 6. 
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Figure 9. - Friction coefficient, displacement thickness, and 
momentum loss thickness in a flat-plate boundary layer. 
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Figure 10. - Velocity profiles in a flat-plate boundary layer. 
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